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AFIT/GAE/AA/ 800-20 


Abstract 

Dynamic  variations  in  pressure  occurring  within  shock  tube  channels 
following  the  exit  of  an  incident  shock  wave  were  investigated.  Three 
channel  sizes  were  investigated.  The  first  size  was  4  in  by  8  in  by  48  in 
(10.16  cm  by  20.32  cm  by  1.22m),  the  second  size  was  4  in  by  1  in  by  25  in 
(10.16  cm  by  2.54  cm  by  63.5  cm),  and  the  third  size  was  1  in  by  1  in  by 
44  in  (2.54  cm  by  2.54  cm  by  1.12  m).  For  the  first  and  second  channels, 
the  shock  tube  diaphragm  pressure  ratio,  ,  ranged  from  6.0  to  20.0. 

The  third  channel  was  used  over  a  range  of  8.0  to  20.0  and  also  to 
determine  the  presence  of  transverse  waves  due  to  an  area  reduction  in 
the  shock  tube.  The  first  channel,  having  the  same  cross-sectional  area 
as  the  shock  tube,  exhibited  pressures  which  approximated  simple  shock 
tube  theory  results.  The  second  and  third  channels  exhibited  complex 
pressure  changes  which  consistently  peaked  at  the  same  time  after  the 
shock  wave  passed  the  pressure  sensors,  independent  of  P^. 
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I.  Introduction 


Under  conditions  of  very  high  entropy,  as  exhibited  in  chemical 
lasers,  undesirable  pressure  waves  can  be  generated  which  cause  the 
efficiency  of  the  system  to  be  degraded.  A  discussion  of  the  phenomenon 
of  high  entropy  generated  pressure  waves  will  provide  the  foundation  for 
this  research  effort. 

A  chemical  laser  utilizes  a  chemical  reaction  to  achieve  the  nec¬ 
essary  atomic  or  molecular  population  inversion  necessary  for  excitation 
energy.  In  the  case  of  a  pulsed  chemical  laser,  these  reactions  can  be 
explosive  in  nature,  generating  a  shock  wave  in  the  gas  medium.  As  the 
shock  wave  propagates  through  the  device,  it  may  reflect  from  a  fixed 
surface  as  a  shock  wave  or  from  a  free  surface  as  an  expansion  fan.  De¬ 
pending  upon  the  duration  and  intensity  of  these  waves,  as  well  as  the 
internal  configuration  of  the  device,  the  pressure  waves  can  generate  a 
very  complex  pressure  distribution  which  directly  affects  the  power  effi¬ 
ciency  and  pulse  rate  of  the  laser. 

Background 

In  an  effort  to  understand  the  phenomenon  of  pressure  waves  in  a 
confined  cavity,  shock  tubes  have  historically  been  used  as  a  shock  wave 
generating  source.  A  previous  approach  in  attempting  to  eliminate  both 
shock  waves  and  expansion  fans  from  propagating  back  into  the  shock  tube 
was  done  by  Jones  and  McCallun/.  They  successfully  investigated  using 
flat  plates  and  grids  as  reflection  eliminators,  causing  the  negative  and 
positive  pressure  waves  to  cancel  each  other  prior  to  reentering  the 
shock  tube. 
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There  may  be  conditions  under  which  it  is  not  possible  to  elimi¬ 
nate  pressure  waves,  either  with  reflection  eliminators  or  by  other  means. 
It  then  becomes  necessary  to  consider  the  impact  of  this  physical  phenom¬ 


enon  on  the  system  of  interest. 

Objective 

The  objective  of  this  research  will  be  to  obtain  qualitative  and 
quantitative  measurements  of  pressure  changes  in  a  channel  following  the 
exit  of  an  incident  shock  wave.  This  is  to  be  accomplished  by  determin¬ 
ing  intensity  and  duration  of  shock  waves  and  expansion  fans  for  the 
defined  conditions. 

Approach 

Tests  will  be  conducted  in  the  Air  Force  Institute  of  Technology 
4  in  by  8  in  by  20  ft  (10.16  cm  by  20.32  cm  by  6.10  m)  shock  tube  facility. 
Ambient  temperature  air  will  be  used  for  both  driven  and  driver  gas.  Shock 

tube  performance  estimates  will  be  made  using  the  basic  relationships  of 

8  9 

Liepman  and  Roshko  and  Shapiro  . 

The  three  test  sections  to  be  used  in  the  experiment  will  be  1  in 

by  1  in,  4  in  by  1  in,  and  4  in  by  8  in  (2.54  cm  by  2.C4  cm,  10.16  cm  by 

2.54  cm,  and  10.16  cm  by  20.32  cm)  in  cross  section.  Surface  mounted 

pressure  transducers,  located  in  each  test  section  at  known  locations, 

will  be  used  to  obtain  the  data. 

Before  discussing  the  results  of  the  investigation,  it  is  necessary 
to  first  consider  a  few  theoretical  aspects  of  shock  tubes  and  to  discuss 
the  unique  features  of  the  shock  tube  and  instrumentation  used  for  this 
research.  Following  these  aspects,  an  analysis  of  the  results  obtained 
from  the  three  test  sections  will  be  presented. 


I I .  Theoretical  Considerations 


To  simplify  the  experimental  conditions  within  the  shock  tube,  the 
following  assumptions  are  made:  (1)  1-D  flow  shock  wave  prior  to  area 
change,  2-D  flow  afterward,  (2)  1-D  flow  expansion  waves,  (3)  air  on  both 
sides  of  diaphragm  is  stationary  and  at  ambient  temperature  prior  to  dia¬ 
phragm  rupture,  and  (4)  air  is  considered  a  perfect  gas  with  y  =  1.4. 

In  order  to  increase  the  effective  Mach  number  range  of  the  shock 
tube  above  that  obtained  by  simply  pressurizing  the  driver  section,  two 
approaches  are  used:  (1)  the  driven  section  is  evacuated,  increasing 
the  diaphragm  pressure  ratio,  P^ ,  and  (2)  the  cross-sectional  area  of 
the  shock  tube  is  reduced  through  the  test  section  region. 

The  operation  of  a  simple  shock  tube  can  best  be  described  with 
the  aid  of  Fig  1.  Prior  to  diaphragm  rupture,  the  driver  section  is  a 
region  of  high  pressure  (P4)  and  the  driven  section  is  a  region  of  low 
pressure  (P^.  At  the  time  of  rupture,  t  =  0,  a  shock  wave  is  initiated 
which  moves  through  the  region  of  low  pressure  while  an  expansion  fan 
moves  into  the  high  pressure  region.  The  contact  surface  between  re¬ 
gions  2  and  3  represents  the  boundary  between  the  gases  that  were 
initially  on  either  side  of  the  diaphragm.  Regions  2  and  3  differ  in 
temperature  but  have  equal  pressures  and  velocities. 

After  the  shock  wave  exits  the  end  of  the  tube,  at  t  =  t g,  an 
expansion  fan  begins  to  propagate  back  into  the  shock  tube  at  sonic  velo¬ 
city,  with  respect  to  the  gas  medium  into  which  it  is  moving.  For  the 
condition  where  the  flow  behind  the  shock  wave  is  equal  to  or  greater 
than  the  speed  of  sound,  Mg  -  1,  the  expansion  fan  will  not  propagate 
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Fig  1,  Simple  Shock  Tube  Dynamic  Processes 


into  the  shock  tube.  With  the  condition  that  air  at  uniform  temperature 
is  used  for  both  driver  and  driven  gas,  the  sonic  flow  condition  behind 
the  shock  wave  is  achieved  when  Mg  =  2.07  ^ . 

The  downstream  pressure,  P-j ,  remains  constant  after  the  exit  of 
the  shock  wave  from  the  shock  tube  if  the  shock  wave  expands  into  an 
infinite  region.  As  the  shock  wave  attempts  to  distribute  its  energy 
over  its  entire  surface,  the  pressure  difference  across  the  shock  dimin¬ 
ishes  as  the  shock  surface  becomes  increasingly  large. 

For  the  shock  tube  with  a  reduced  area  test  section,  the  dynamic 
processes  behave  in  a  more  complex  manner,  as  shown  in  Fig  2.  As  the 
incident  shock  wave  moves  into  the  reduced  area  test  region.  Fig  2b,  a 
process  occurs  which  is  described  by  Tambra  10  as  "coalescing"  trans¬ 
verse  waves.  The  area  change  causes  the  formation  of  these  waves  which 
move  transverse  to  the  initial  flow,  yet  also  remain  attached  to  the 
incident  shock  wave  as  it  moves  down  the  test  section.  As  the  transverse 
waves  move  across  the  incident  shock  surface  they  tend  to  strengthen  the 
incident  shock  but  diminish  in  strength  over  distance.  Along  with  the 
development  of  the  transverse  waves  as  a  result  of  the  area  reduction, 
part  of  the  incident  shock  wave  is  also  reflected  off  the  ramps  and  prop¬ 
agates  back  upstream. 

Following  the  exit  of  the  incident  shock  wave  from  the  reduced 
area.  Fig  2c, the  shock  wave  will  continue  to  move  downstream  as  an  expan¬ 
sion  fan  moves  upstream  into  the  reduced  area  region.  If  the  expansion 
into  the  larger  part  of  the  test  section  does  not  dissipate  the  shock  wave, 
pressure  P -j  will  remain  higher  than  P^.  This  condition  reduces  the 
strength  of  the  expansion  fan  that  moves  back  into  the  reduced  area  region. 
When  the  expansion  fan  reaches  the  upstream  end  of  the  reduced  area  test 
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section,  it  reflects  off  the  "free  surface"  of  the  larger  area  as  a  series 
of  compression  waves.  Fig  2d.  These  compression  waves  combine  as  they 
propagate  through  the  test  section  to  become  a  stronger  single  shock  wave. 
Finally,  as  the  incident  shock  wave  leaves  the  end  of  the  shock  tube,  another 
expansion  fan  moves  upstream  at  a  velocity  determined  by  the  ratio. 

The  theoretical  method  used  to  determine  enhanced  shock  Mach  numbers 
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after  an  area  change  was  to  apply  Whitham's  rule  *  .  This  rule  is  an  empir¬ 

ical  relationship  which  assumes  instantaneous  Mach  number  changes  following 
abrupt  area  changes.  The  method  is  elaborated  in  Appendix  B.  For  the  series 
of  test  runs,  a  computer  program  was  written  to  calculate  expected  Mach 
numbers  for  the  area  reductions  considered,  given  incremental  ratios. 
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III. 


Shock  Tube  System 

The  shock  tube  used  in  the  research  project  was  the  Air  Force 
Institute  of  Technology  4  in  by  8  in  by  20  ft  (10.16  cm  by  20.32  cm  by 
6.10  m)  shock  tube"*’  The  driver  section  is  4  ft  (1.22  m)  in  length. 

This  section  is  mounted  on  a  sliding  and  traversing  mechanism  which  per¬ 
mits  access  to  both  sections  for  cleaning.  The  driven  section  is  16  ft 
(4.88  m)  in  length,  constructed  of  three  straight  4  ft  (1.22  m)  sections 
and  a  changeable  4  ft  (1.22  m)  test  section  in  which  the  internal  dimen¬ 
sions  can  be  altered.  A  shock  dump  tank  which  can  be  attached  to  the  end 
of  the  shock  tube  will  be  described  later. 

Gas  pressure  in  the  driver  section  is  monitored  by  a  calibrated  0-200 

O 

in  Hg  (0  -  677.3  kN  -  m  )  bourdon  gage.  Vacuum  in  the  driven  section  is 
monitored  by  a  30  in  Hg  (101.59  kN  -  m  )  manometer  tube.  Air  was  used 
as  the  driver  and  driven  gas.  The  air  supply  was  the  filtered,  dry  com¬ 
pressed  air  from  the  building  compressors.  Maximum  air  supply  pressure 

O 

was  approximately  228  in  Hg  (772.1  kN  -  m  )  gage.  Vacuum  in  the  driven 
section  was  achieved  by  a  Heraeus  Model  E  70  Vacuum  pump.  Maximum  attain¬ 
able  vacuum  with  the  dump  tank  attached  was  approximately  22  in  Hg  (74.5 
kN  -  m  ).  Figure  3  is  a  diagram  of  the  shock  tube  system. 

The  diaphragm  used  to  maintain  the  pressure  differential  in  the 
shock  tube  was  Mylar,  used  in  either  single  or  multiple  thicknesses. 

Previous  reports  indicated  using  no  more  than  four  sheets  of  Mylar  at  a 
3  11 

time  ’  .  In  this  investigation,  no  difficulties  were  experienced  using 

5  or  6  sheets  of  Mylar,  1x10"^  or  2x10"^  in  (2.54x10"^  or  5.08xl0"5m)  in 
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thickness. 

Myl ar  is  much  more  resistant  to  shattering  upon  rupture  than  other 

non-metallic  materials  and  hence  provides  more  consistent  shock  tube 

performance.  Furthermore,  it  is  necessary  to  establish  the  pressure 

differential  as  near  the  bursting  strength  of  the  diaphragm  as  possible 

to  achieve  the  uniform  shock  velocity  downstream  in  the  shortest  possible 
4  5 

distance  ’  .  The  Mylar  thickness  required  to  maintain  the  pressure 

-5  -7  -2 

differential  during  this  experiment  was  4x10  in/in  Hg  (3x10  m/kN  -  m  ). 
Previous  investigations  have  recommended  greater  thicknesses  to  maintain 
pressure  differentials  with  this  shock  tube.  The  diaphragm  was  ruptured 
on  command  by  a  remotely  controlled  pyramidal  probe. 

Test  Sections 

Three  4  ft  (1.22  m)  test  sections  were  used  during  the  research. 
Significant  dimensions  for  the  test  sections  are  shown  in  Fig  4.  Test 
section  1  was  used  as  a  straight,  unchanged  section.  Sections  2  and  3 
were  modified  with  pairs  of  rectangular  wood  blocks  which  provided  a  1  in 
(2.54  cm)  high  open  channel  in  the  middle  of  the  test  sections.  The 
straight  taper  ramps  a  and  b,  also  in  pairs,  were  positioned  at  the 
driver  end  of  the  appropriate  test  section.  For  test  section  2,  ramp  a 
was  used  to  reduce  the  height  from  8  in  (20.32  cm)  to  1  in  (2.54  cm). 

Test  section  3  also  used  ramp  a  to  reduce  the  height  from  8  in  (20.32  cm) 
to  1  in  (2.54  cm)  and  used  ramp  b  to  reduce  the  width  from  4  in  (10.16  cm) 
to  1  in  (2.54  cm). 
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Upon  initiation  of  the  shock  in  a  constant  area  shock  tube  with 
a  closed  end  on  the  driven  section,  the  shock  propagates  to  the  end  of 
the  driven  section  and  is  reflected  back  through  the  test  region  as  a 
normal  shock.  In  order  to  provide  the  condition  of  a  shock  wave  exiting 
a  confined  channel  and  initiating  the  expansion  fan  process,  the  end  of 
the  shock  tube  was  equipped  with  a  shock  dump  tank.  The  tank  is  cylind¬ 
rical  in  shape,  32  in  (81.28  cm)  in  diameter  and  6  ft  (1.83  m)  in  length. 

The  tank  provides  the  additional  benefits  of  allowing  the  driven  section 
to  be  evacuated,  dissipating  the  shock  waves  from  the  shock  tube  and 
delaying  any  shock  reflections  which  might  propagate  into  the  test  sec¬ 
tion.  The  dump  tank  was  designed  in  a  modular  fashion  to  permit  attach¬ 
ment  to  the  end  of  all  three  test  sections. 

Electronic  Instrumentation 

The  electronic  instrumentation  used  for  this  research  included 
pressure  transducers,  amplifiers,  oscilloscopes,  digital  timers,  a  fre¬ 
quency  oscillator,  and  a  time  delay  computer.  The  electronics  arrange¬ 
ment  is  presented  diagrammatical^  in  Fig  5.  Sensors  A  and  B,  Endevco 
8501-500  pressure  transducers,  were  used  to  detect  the  incident  shock 
speed  and  trigger  the  time  delay  computer.  Sensors  C  and  D,  Endevco  8506- 
15  and  8506-50  pressure  transducers,  respectively,  were  used  to  measure 
the  pressure  fluctuations  in  the  test  sections  caused  by  the  shock  waves  and 
expansion  fans.  All  transducers  were  surface  mounted  in  the  upper  side  of 
the  test  sections.  For  the  three  test  configurations,  the  appropriate 
distances  between  transducers  are  given  in  Table  1. 
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Fl9  5'  Electr°nic  Arrangement 


The  time  delay  computer  developed  by  Weber11  was  used  to  trigger 
the  sweep  of  the  oscilloscopes.  This  was  accomplished  with  the  use  of 
an  internal  counter  used  to  count  the  number  of  pulses  generated  by  a 
10  Mhz  ( f i )  oscillator,  starting  from  zero,  as  the  incident  shock  wave 
traveled  the  distance  d-j  from  sensor  A  to  sensor  B. 

The  same  number  of  pulses  was  then  counted  down  to  zero  at  the  f2 
frequency  rate  as  the  shock  moved  a  distance  d2>  At  that  time  an  output 
signal  was  generated  which  externally  triggered  the  oscilloscopes.  For 
a  constant  shock  speed,  the  frequency  ratio  f-|/f2  is  equal  to  the  dis¬ 
tance  ratio  d2/d-|  and  is  independent  of  shock  Mach  number.  The  f2 
oscillator  was  therefore  adjusted  to  initiate  the  sweep  of  the  oscillo¬ 
scopes  with  the  shock  at  a  position  in  the  test  section  just  prior  to 
arrival  at  sensor  C. 
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IV.  Results 


General 

The  emphasis  of  this  research  effort  was  to  obtain  pressure  measure¬ 
ments,  following  the  passage  of  an  incident  shock  wave,  in  three  different 
shock  tube  test  sections,  differing  in  cross-sectional  area.  In  addition, 
test  section  3  was  used  to  verify  the  presence  of  "coalescing"  tranverse 
waves  following  an  area  reduction  in  the  shock  tube,  as  described  by 
Tambra^. 

In  Tambra's  work,  he  investigated  the  change  in  shock  wave  velocity 
due  to  the  transverse  waves  following  abrupt  area  changes  over  a  P41 
pressure  ratio  of  4.0  to  7.0.  He  observed  that  as  the  shock  travelled 
through  the  reduced  area  section,  it  varied  in  speed,  usually  increasing 
to  a  maximum  velocity  and  then  decreasing.  The  distance  downstream  of 
the  area  change  at  which  the  maximum  velocity  occurred  changed  very  little 
with  changes  in  P41  but  was  dependent  upon  the  area  ratio  of  the  reduced 
area  to  the  initial  area.  As  the  area  ratio  decreased,  the  maximum  speed 
location  point  moved  toward  the  point  of  area  change. 

Toepler  schlieren  photography®  was  used  to  determine  the  presence  of 
transverse  waves  immediately  behind  the  incident  shock  in  test  section  3. 
Although  the  transverse  waves  were  present  at  P^  ratios  similar  to 
Tambra's  study,  they  were  weaker.  This  weakness  may  be  attributable  to  the 
long  distance  downstream  of  the  area  change  that  the  photographs  were 
taken  and  to  the  tapered  ramps.  The  long  distance  involved  allowed  the 
waves  to  diminish  in  strength  whereas  the  ramps  may  have  provided  a 
different  effect  on  the  strength  of  the  shock  in  the  reduced  area  section 
than  occurred  for  the  abrupt  area  change  investigated  by  Tambra. 


Measured  Mach  numbers  corresponded  very  well  with  results  predicted 
by  simple  shock  tube  theory  and  Whitham's  rule  for  the  three  test  sections. 
For  the  determination  of  shock  Mach  number,  M  ,  over  the  range  investi¬ 
gated,  Whitham's  rule,  used  for  test  sections  2  and  3,  is  almost  as  accur¬ 
ate  as  the  theory  of  straight  shock  flow  used  for  test  section  1.  Results 
for  the  three  sections  are  given  in  Table  2.  With  the  use  of  different 
thicknesses  of  Mylar  sheets  to  produce  the  same  total  thickness  required 
for  a  given  P^ ,  Mg  varied  less  than  0.6%  over  the  P^  range  investigated. 

Figure  6  shows  that  the  initial  transient  response  of  sensor  C  to  the 
passage  of  the  incident  shock  wave  was  not  always  the  same  with  the  same 
test  conditions.  One  trace  indicates  a  much  faster  and  larger  initial 
response  to  the  shock  wave  than  the  other  trace.  However,  the  traces  have 
the  same  value  after  approximately  50  usee  following  the  shock  wave  pass¬ 
age,  indicating  the  same  test  section  pressure. 

For  test  sections  1  and  2,  the  P^  ratios  of  6.0  and  7.0  were  the  only 
runs  in  which  the  driven  section  was  open  to  the  atmosphere  and  the  dump 
tank  was  not  attached.  The  other  tests  were  conducted  with  P^  less  than 
atmospheric  pressure.  In  order  to  compare  test  results  of  the  different 
P-j  pressures, transducer  output  signals  were  nondimensional ized  by  the 
initial  driven  pressure  used  for  each  run.  This  nondimensional ized  test 

I  I 

section  pressure  is  ?2  for  test  section  1  and  Pg  for  test  sections  2  and 

I  I 

3.  The  ratios  Pg  and  Pg  are  equal  to  P and  Pg-j ,  respectively,  for  the 
condition  when  the  incident  shock  wave  passes  the  sensor. 

Test  Section  1 

I 

Figures  7  and  8  show  the  Pg  responses  recorded  by  sensors  C  and  D 

I 

in  test  section  1.  As  P^  increases,  the  initial  equilibrium  Pg  ratio 


Table  2.  Expected  vs  Measured  Results  for  M$,  and  Pg-j 


Test  Section 

P41 

P-j(psia) 

M  M  % 

se  sm 

error 

P21e 

P21m 

%  error 

1 

6.0 

14.38 

1.46  1.39 

5.0 

2.320 

2.087 

12 

7.0 

14.21 

1.50  1.43 

4.9 

2.458 

2.219 

11 

8.0 

11.31 

1.54  1.48 

4.1 

2.600 

2.389 

8.8 

10.0 

9.35 

1.61  1.54 

4.5 

2.857 

2.600 

9.9 

15.0 

6.40 

1.74  1.66 

4.8 

3.366 

3.048 

10 

20.0 

5.41 

1.83  ' 

1.76 

4.0 

3.740 

3.447 

8.5 

Test  Section 

P41 

P^psia) 

Mse 

M 

sm 

%  error  P,, 
ole 

P61m 

%  error 

2 

6.0 

14.23 

1.60 

1.64 

-2.4 

2.820 

2.971 

-5.1 

7.0 

14.23 

1.66 

1.68 

-1.2 

3.048 

3.126 

-2.5 

8.0 

10.67 

1.71 

1.76 

-2.8 

3.245 

3.447 

-5.9 

10.0 

9.35 

1.79 

1.84 

-2.7 

3.571 

3.783 

-5.6 

14.0 

6.40 

1.93 

2.00 

-3.5 

4.179 

4.500 

-7.1 

20.0 

4.43 

2.07 

2.21 

-6.3 

4.832 

5.531 

-12 

Test  Section  P^ 

P-j(psia) 

Mse 

^sm 

%  error 

P61e 

P61m 

%  error 

3  8.0 

11.42 

1.73 

1.84 

-6.0 

3.325 

3.783 

-12 

10.0 

9.45 

1.82 

1.97 

-7.6 

3.698 

4.361 

-15 

14.0 

6.58 

1.95 

2.18 

-11 

4.270 

5.378 

-21 

20.0 

5.60 

2.09 

2.43 

-14 

4.929 

6.722 

-27 

3 


also  increases.  However,  in  comparing  with  P21m  from  Table  2,  P2  in¬ 
creases  at  a  faster  rate  than  P^m.  P2  a^so  maintains  a  constant  value 
for  a  longer  period  of  time  and  does  not  exhibit  as  rapid  a  decrease  as 
P^  increases.  The  shock  wave  pressure  ratio,  P^,  and  the  Mach  number 
of  the  flow  behind  the  shock  wave,  M2,  increase  with  increasing  P^. 

These  conditions  slow  the  movement  of  the  expansion  fan  into  the  test 
section  and  cause  a  delay  in  pressure  decrease,  with  sensor  C  exhibiting 
the  expansion  fan  effect  later  than  sensor  D.  Table  3  compares  the 
theoretical  arrival  times  of  the  1-D  expansion  fan  phenomenon  to  the 
actual  arrival  times  for  the  different  P^  ratios  in  test  section  1. 

The  times  determined  theoretically  assume  that  when  the  constant  velocity 
incident  shock  wave  travels  to  the  end  of  the  reduced  area  section,  the 
change  in  boundary  conditions  causes  an  expansion  fan  to  immediately 
begin  to  propagate  back  into  the  reduced  area  section.  The  times  in¬ 
clude  shock  wave  travel  from  the  sensors  to  the  end  of  the  area  reduction 
and  the  expansion  fan  travel  back  to  the  sensors.  Due  to  the  differences 
in  velocities  of  the  shock  wave  and  the  contact  surface,  and  the  length  of 
the  shock  tube,  the  contact  surface  was  determined  not  to  arrive  at  the 
sensor  locations  until  after  the  expansion  fan  process  had  occurred.  The 
observation  that  the  actual  travel  times  are  usually  longer  than  the 
theoretical  travel  times  indicates  that  the  expansion  fan  probably  does 
not  start  to  propagate  immediately  after  the  shock  wave  exits  the  channel. 

I 

Rather  than  P2  exhibiting  a  continuous  gradual  decrease  as  a  result  of 
the  expansion  fan  process,  P2'  decreases  to  an  intermediate  level  before 
decreasing  to  its  final  value.  This  process  is  more  apparent  with 
sensor  D.  This  suggests  that  the  initial  expansion  fan  that  enters  the 
test  section  establishes  an  equilibrium  state  which  is  subsequently 

followed  by  an  additional  expansion  fan. 
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Table  3.  Theoretical  vs  Experimental  Propagation  Times 
for  Expansion  Fan,  Test  Section  1 _ 

P,,  Actual  Arrival  Time  (msec)  Theoretical  Arrival  Time 

(msec) 


C 

D 

C 

D 

6.0 

1700 

600 

1672 

557 

7.0 

1700 

600 

1738 

579 

8.0 

1900 

760 

1851 

617 

10.0 

2200 

750 

2025 

675 

15.0 

2700 

1150 

2543 

848 

20.0 

(obscured) 

1500 

3249 

1083 

Test  Section  2 

I 

Figures  9  and  10  are  the  Pg  responses  for  sensors  C  and  D  in 
test  section  2.  These  curves  show  the  same  trend  as  P2  in  test  section  1 

I 

for  Pg  to  increase  faster  than  Pglm  as  P^-j  increases.  Immediately  after 

I 

the  passage  of  the  shock  wave  across  the  sensors,  a  decrease  in  Pg  is 
followed  by  an  increase  and  decrease.  For  both  sensors  C  and  D,  this 
pressure  spike  appears  to  be  related  to  P41  because  the  Pg  increase  be¬ 
comes  larger  as  P41  also  increases.  With  sensor  C,  the  increase  also 
occurs  at  a  progressively  later  time,  moving  from  .45  msec  to  1.16  msec 
after  the  incident  shock  wave,  over  the  range  of  P^.  For  sensor  D,  the 

I 

occurrence  of  the  Pg  increase  changes  minimally. 

For  the  lower  values  of  P41 ,  the  curves  are  similar  to  test  section 
1  except  for  the  initial  pressure  disturbance.  For  P^  of  6.0  and  7.0, 
the  intermediate  pressure  plateau  is  present  and  is  also  present  in  the 
curves  for  P^  of  8.0  and  10.0,  but  with  the  latter  two  it  is  followed 
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by  a  Pg  increase  instead  of  a  decrease.  For  the  range  of  P^  from  8.0 
to  20.0,  the  pressure  rise  continues  to  occur  with  its  maximum  value 
achieved  at  6  msec  following  the  passage  of  the  shock  wave.  This  press¬ 
ure  rise  begins  sooner  after  the  passage  of  the  shock  wave  as  P^  in¬ 
creases. 

The  Pg  rise  at  6  msec  coincides  with  the  time  it  would  theoretically 
take  for  an  expansion  fan  to  propagate  to  the  front  of  the  reduced  area 
section  and  return  as  a  shock  wave  for  ratios  of  8.0  and  10.0,  assum¬ 
ing  a  downstream  pressure  of  P.] .  This  observation  tends  to  support  the 
idea  that  for  these  conditions,  P -j  is  approximately  equal  to  P^.  The 
gradual  change  from  pressure  decrease  to  pressure  increase  over  the  range 
of  P^i  is  not  understood,  but  for  P^  =  20,  an  expansion  fan  cannot  enter 
the  test  section  because  I^l.  Were  the  pressure  rise  a  result  of  the 
expansion  fan  becoming  a  shock  wave  at  the  upstream  end  of  the  test  section, 
the  longer  time  required  to  propagate  upstream  would  delay  its  occurrence 
and  the  wave  would  not  be  present  at  all  for  P41  =  20.  The  very  fact  this 
is  not  what  happens  indicates  the  presence  of  a  very  complex  system  of 
pressure  changes  occurring  in  the  shock  tube  with  pressure  increases  be¬ 
coming  more  dominant  for  increases  in  P^,  The  complexity  of  pressure 
changes  in  the  test  section  is  apparently  due  to  the  reduction  of  cross- 
sectional  area  in  the  test  section. 

Test  Section  3 

I 

Figures  11  and  12  are  the  Pg  responses  of  sensors  C  and  D  in  test 
section  3.  From  Table  2  there  is  a  more  significant  difference  between 
measured  and  expected  values  of  Pg-j  in  section  3  than  in  section  1  or  2. 

I 

Despite  this  larger  difference,  the  initial  value  of  Pg  in  section  3 
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corresponds  more  closely  to  the  measured  value  of  Pg-j  over  the  range  of 
P^I  than  in  section  1  or  2.  As  with  section  2,  a  pressure  disturbance 
occurs  shortly  after  the  passage  of  the  incident  shock  wave.  However, 
the  time  the  disturbance  occurs  does  not  change  as  P^  increases.  For 
P^-j  ratios  of  8.0  and  10.0,  the  pressure  curves  are  very  similar  to 
those  of  section  2.  Even  though  Pg-|  is  larger  in  section  3  than  sec- 

I 

tion  2  due  to  the  greater  area  reduction,  the  maximum  value  of  Pg  in 
section  3  is  significantly  less  than  in  section  2.  This  suggests  that 
as  the  distance  changes  downstream  from  the  area  change,  the  peak  press¬ 
ure  achieved  at  that  point  in  the  test  section  will  also  change.  The 

I 

peak  value  for  Pg  still  occurs  at  approximately  6  msec  after  shock  wave 
passage. 
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Sensor  C,  Test  Section  1  traces:  P-,=6.0,  50  ysec/div 
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vs  Time.  Test  Section  1,  Sensor  C 


1  2  3  4  5  6  P4l“ 

Time  (msec) 


Fig.  10  Pg1  vs  Time,  Test  Section  2,  Sensor  D 
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V.  Conclusions 


The  results  of  this  investigation  indicate  pressure  variations  in 
a  channel,  after  the  exit  of  an  incident  shock  wave,  are  significantly 
influenced  by  changes  in  length  and  cross-sectional  area.  In  a  constant 
area  shock  tube,  for  a  P^  range  of  6.0  to  20.0,  the  expansion  fan  pro¬ 
cess  is  the  dominant  pressure  influence.  In  a  shock  tube  with  an  area 
change,  the  pressure  variation  becomes  a  much  more  complex  sequence  of 
pressure  changes  following  the  incident  shock  exit. 

In  test  section  1,  the  expansion  fan  process  following  the  exit  of 
the  incident  shock  wave  does  not  appear  to  be  a  continuous  process,  but 
one  that  occurs  in  several  quasi-continuous  increments.  In  addition,  the 

I 

initial  P21  ratio,  as  determined  from  P 2  ,  becomes  increasingly  larger 
than  expected  as  P41  increases. 

In  test  sections  2  and  3,  the  expansion  fan  process  is  observable 
at  the  lower  P41  ratios  investigated  and  behaves  in  a  similar  manner  to  the 
process  in  test  section  1.  As  P^  increases,  large  pressure  increases 
become  the  dominant  pressure  variation  within  the  test  section.  The 
pressure  increases  consistently  occur  at  approximately  the  same  time  after 
the  passage  of  the  incident  shock  wave  regardless  of  P^.  The  intensity 
of  the  pressure  increases  is  a  function  of  length  and  cross-sectional 
area  of  the  test  sections.  In  test  section  2,  as  in  test  section  1,  the 

I 

initial  Pg^  ratio,  determined  from  Pg  ,  becomes  increasingly  larger  than 
expected  as  P^  increases.  Due  to  the  complex  nature  of  the  flow  behind 
the  incident  shock  wave  in  a  shock  tube  with  a  reduced  area  test  section, 
additional  investigation  is  necessary  to  explain  these  occurrences. 
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VI.  Recommendations 


There  are  several  recommendations  that  can  be  proposed  for  addi¬ 
tional  research  in  this  area  of  study: 

1.  Use  4  in  (10.16  cm)  wide  test  section  for  optical  investigation 
of  transverse  wave  phenomenon  in  conjunction  with  pressure  sensors  posi¬ 
tioned  along  the  entire  length  of  the  test  section.  This  would  permit 
determination  if  the  pressure  fluctuations  may  be  correlated  with  the 
transverse  wave  phenomenon. 

2.  Investigate  study  of  pressure  fluctuations  in  reduced  area 
sections  which  extend  through  the  entire  test  section.  This  study  would 
eliminate  any  effects  caused  by  a  change  in  boundary  conditions  within 
the  test  section  prior  to  the  exit  of  the  shock  wave  from  the  end  of  the 
test  section. 

3.  Investigate  and  compare  results  in  similar  test  sections  with  a 
shock  wave  generated  by  a  chemical  reaction  or  by  electric  discharge. 
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Appendix  A 


Experimental  Procedure 

Initial  configuration: 

-  P.j  and  gage  shutoff  valves  -  open 

-  P^  pressurization  and  P^  bleed  valves  -  closed 

-  Plunger  and  plunger  bleed  valves  -  closed 

-  All  nonessential  valves  -  closed 

1.  Clean  shock  tube  and  install  appropriate  thicknesses  of  diaphragm 
material . 

2.  Zero  all  digital  timers  and  oscilloscopes  and  set  sweep  trigger  on 
oscilloscopes. 

3.  Adjust  frequency  for  desired  shock  location. 

4.  Open  vacuum  solenoid  valve,  start  vacuum  pump  and  open  plunger  valve. 

5.  Open  P^  pressurization  valve. 

6.  When  desired  diaphragm  pressure  ratio  is  reached,  close  vacuum  sole¬ 
noid  valve. 

7.  Close  P^  and  P^  gage  shutoff  valves  and  remotely  rupture  diaphragm. 

8.  Shut  off  vacuum  pump,  close  plunger  valve,  open  plunger  bleed  valve 
and  P^  gage  shutoff  valve  and  open  P^  pressurization  valve  more  to  allow 
internal  pressure  to  reach  atmospheric  level. 

9.  Record  digital  timers  At^  and  Atg  .  Photograph  oscilloscope  traces. 

10.  Mark  photographs  to  correspond  with  numbered  test  run. 

11.  When  P.|  gage  pressure  reaches  atmospheric  pressure  close  P4  pressuri¬ 
zation  valve  and  plunger  bleed  valve  and  open  P^  gage  shutoff  valve. 

12.  Open  shock  tube,  reset  plunger  probe  and  clean  out  shock  tube  debris 


as  necessary. 


Appendix  B 


Frequency:  adjustable 

Freq  Multiplier  Switch:  x  1M,  x  1 OK 

Symmetry:  Normal 

Gen  Mode:  Continuous 

Signal  Output:  Pulse  Out  (TTL) 


Hewlett-Packard  Model  5325B  Universal  Counter 
(At-|  configuration) 

Time  Base:  .1  ysec 
Selection:  T.I.  A  to  B 
Channel  A:  Atten:  x  10 
AC 

Slope:  (  -  ) 

Channel  B:  Atten:  x  10 
AC 
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Slope:  (  +  ) 

(At2  configuration) 

Time  Base:  .1  ysec 
Selection:  T.  I.  A  to  B 
Channel  A:  Atten:  x  10 
AC 

Slope:  {  -  ) 

Channel  B:  Atten:  x  10 
DC 

Slope:  (  -  ) 

Tektronix  Model  434  Storage  Oscilloscopes 
Mode:  single 
Trigger  Source:  external 
Slope:  (  +  ) 

Level :  (  +  ) 

Ext  Attenuation:  1:1 
LF  Rej:  off 
HF  Rej :  off 
Coupling:  DC 

Electro  Instruments  Model  A20B-1  Amplifier  Specifications 
Full  Scale  Voltage:  +  10  vdc  or  peak  AC  to  20  Khz 
Linearity:  DC:  +  .01%  max  (100  ohm  min  load) 

AC:  DC  slew  shift  not  more  than  +  .01%  of  peak  AC 
up  to  1  Khz 

Freq  Response:  +  1%  DC  to  2  Khz.  1  dp  to  20  Khz, 

3  db  to  50  Khz 
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APPENDIX  C 
Data  Reduction 


Whitham's  Rule 


where 


2 

K(M)  *  2{(1  +  — —  H-  )  (2y  +  1  +  M-2))"1 

y  +  1  y 


For  A-|  =  32  sq  In  (206.5  sq  cm),  A?  =  4  sq  in  (25.8  sq  cm)  and 
M$  *  1.60: 


y2  *  .44681 
K(M)  *  .43398 


Therefore  6M  =  .185 

and  Ms  *  1.60  +  .185  *=  1.785  following  the 
area  change. 
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Pressure  Ratio  Calculatl on 


P.  *  p/i  x  c  +  pa  “  P_  x  1.002  +  Pa 
4  g  a  g  a 


P1  "  Pa-hl.r 


For  Pa  =  29.02  In  Hg  (98.277  kN  -  m"£), 

a 

Pg  *  145  In  Hg  (491.047  kN  -  m"2),  and 
h1  r  «  7.3  In  Hg  (24.722  kN  -  m-2): 


P4  =  145  x  1.002  +  29.02  =  174.3  In  Hg  (590.272  kN  -  nfZ) 

P1  «  29.02  -  7.3  =  21.72  In  Hg  (73.555  kN  -  nf2) 

n  —  174.3  _  on 
P4i  rnn  8*° 


Fundamental  Pressure  Ratio 


P41  *  P21  fl  - 


Clf4  -  I)  <Va4)  (p21  -  ')  -2y4/  (y4  -  1) 


r,  -  iy, 


Shock  Mach  Number  (Mg) 


Y1  "  1  Yi  +  1  5 

\  *  {-V  +  TyT1’?!1 


Mach  Number  of  flow  behind  shock  (Mg) 


1  Yi  +  1  Yi  -  1  -.5 

—  (p-.  -  n  (p_.  (  -U —  +  -U —  p_.n 


M2  *  ”y7  (P21  ‘  1 ^  {P21  * 


P21^ 


*  s 


Normalized 


P2  (P2*) 

P-  *  P  +  P„  where  P„  is  determined  by  vertical  displace- 
c  &  g  g 

ments  of  oscilloscope  traces  from  the  0  dc  line. 

V  *  VP1 


For  P. 


29.02  In  Hg  ■  14.26  psia  (98.277  kN  -  m'2)  and 

5  dlv  x  100  mv/div  *  25.45  mv/psi 
19.65  psig  (135.42  kN  -  m"2) 

9.02  in  Hg  •  4.43  psia  (30.547  kN  -  m"2) 


14.26  +  19.65 

VZ  ‘  03 


7.65 
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